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We review recent progress in the determination of the subsaturation density behavior of 
the nuclear symmetry energy from heavy-ion collisions as well as the theoretical progress 
in probing the high density behavior of the symmetry energy in heavy-ion reactions 
induced by high energy radioactive beams. We further discuss the implications of these 
results for the nuclear effective interactions and the neutron skin thickness of heavy 
nuclei. 



1. Introduction 

One of the most important properties of a nuclear matter is its equation of state 
(EOS). For symmetric nuclear matter with equal numbers of protons and neutrons, 
its EOS has been relatively well determined from the study of the nuclear giant 
monopolc resonances (GMR)L!Jas well as the measurements of collective flows'^ and 
subthreshold kaon production 1^ in nucleus-nucleus collisions. On the other hand, 
our knowledge on the EOS of the isospin asymmetric nuclear matter with unequal 
numbers of protons and neutrons, especially the part related to the nuclear symme- 
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try energy, is still largely uncertaini I l ^ l " l ' l ° l ^ l " i Although the nuclear symmetry 

energy at normal nuclear matter d ensity is known to be around 30 McV from the 

empirical liquid-drop mass formula^ I I, its valu es at other densities, particularly 

at supra-normal densities, are poorly known i-^^. Advances in radioactive nuclear 

beam facilities provide, however, the possibility to pin down the density dcpcn- 

dence of the nuclear symmetry energy in heavy ion collisions induced by these nuclei 

I4l5l9ll3ll4ll5ll6ll7ll8ll9l20l2'il22|-^3lkl25l26l27l28l29l30l3ll32l33l34l35l36l37l38l39l40l41l4^l43l 

In the present talk, we review recent progress in extracting the information on 
the subsaturation density behavior of the nuclear symmetry ene rgy from heavy-ion 
collisions, especially from the analysis of isospin diffusion data I 44 | 4o | 4d | p^j-ther- 
more, we discuss the implications the constrained symmetry energy has for the 
nuclear effective interactions and the neutron skin thickness of heavy nuclei. In ad- 
dition, we also review the theoretical progress on probing the high density behavior 
of the symmetry energy in heavy-ion reactions induced by high energy radioactive 
beams. 

2. The nuclear symmetry energy 

The EOS of isospin asymmetric nuclear matter, given by its binding energy per 
nucleon, can be generally written as 

E{p, a) = E{p, a = 0) + Esym{p)a^ + 0(«^), (1) 

where p = pn + Pp is the baryon density with p„ and pp denoting the neutron and 
proton densities, respectively; a = {pn — Pp)/{pp + Pn) is the isospin asymmetry; 
E[p, a = 0) is the binding energy per nucleon in symmetric nuclear matter, and 



1 9^g(p, a), 

2 a^ 



Esyraip) - 7^ -7^ |q=0 (2) 



is the nuclear symmetry energy. Neglecting contributions from higher-order terms in 
Eq. Il]) leads to the well-known empirical parabolic law for the EOS of asymmetric 
nuclear matter. For nuclear matter at moderate densities, this parabolic law has 
been verified in results from all many-body theories. As a good approximation, the 
symmetry energy can thus be extracted from the binding energy difference between 
pure neutron matter and symmetric nuclear matter, i.e., Egymip) ~ E{p, a = 1) — 
E{p,a = 0). It should be mentioned that possible presence of higher-order terms 
in a at supra-normal densities could significantly modify the proton fraction in 
/3-cquilibrium neutron-star matter and thus the critical densi ty for the direct Urea 
process that is responsible for faster cooling of neutron stars I ' I "I 

Around the nuclear matter saturation density po-, the nuclear symmetry energy 
Esymip) can be expanded to second-order in density as 

E.,-„(.)^ £,,.(.„) +f(^)+^(^)^ (3) 
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where L and -ftTsym are the slope and curvature parameters of the nuclear symmetry 
energy at po, i.e., 



Q \p=Poi J^sym — ypo Q2g \p=Po- 



^ = '^^0 ^7 I P=P0 ' ^^ sym = 9Po HaT I P=Po ' \V 



The i and Ksym characterize the density dependence of the nuclear symmetry 
energy around normal nuclear matter density, and thus carry important information 
on the properties of nuclear symmetry energy at both high and low densities. In 
particular, the slope parameter L has been found to correlate linearly with the 
neutron-skin thickness of heavy nuclei and thus can in pr inciple be determined 
from measured thickness of the neutron skin of such nuclei l'l4SJ|i3U | oi | oz 06 j4ii30 | 
Unfortunately, because of the large uncertainties in the experimental measurements, 
no reliable information on the slope parameter of the nuclear symmetry energy has 
so far been obtained in this way. 

At the saturation density and around a = 0, the isobaric incompres sibility of 
asymmetric nuclear matter can also be expressed to second-order in a asESE3 

K{a)^Ko + K^,ya'^, (5) 

where Kq is the incompressibility of symmetric nuclear mat ter at the saturation 
density and the isospin-dependent part A'asy ~ ^^sym ~ 6L LLqI characterizes the 
density dependence of the nuclear symmetry energy. Information on Xasy can in 
principle be extracted experimentally by measuring the GMR in neutron-rich nuclei. 

3. Constraining the subsaturation density behavior of the 
symmetry energy from isospin diffusion data 

The symmetry ener gy i s known to affect the isospin diffusion in intermediate-energy 
heavy ion collisions ^^^. Experimentally, the degree of isospin diffusion between the 
projectile nucleus A and the target nucleus B can be studied via ^^ ■ '^"1 

2 vA+B _ V''4+/l _ vS+s 

where X is any isospin-sensitive observable. By construction, the value of Ri is 
1 (—1) for symmetric A-\- A {B + B) reaction. For asymmetric reactions, the value 
of Ri decreases from 1 to about zero as the colliding nuclei approaches isospin 
equilibrium. In the NSCL/MSU experiments with A = ^^"^Sn and B = ^^^Sn at a 
beam energy of 50 MeV/nucleon and an impact parameter about 6 fm, t he isospin 
asymmetry of the projectile-like residue was used as the isospin tracer Xl^. Using 
an isospin- and momentum-dependent IBUU04 transport model with experimental 
free-space or in-medium nucleon-nucleon (NN) cross sections, the depende nce of 
Ri on the nuclear symmetry energy in these reactions was studied in Refs. Ill 
for different isospin- and momentum-dependent MDI interactions '^^. Shown in 
Fig. [1] are the results for the degree of the isospin diffusion 1 — Ri as functions of 
iCasy (left panel) and L (right panel). The shaded band in Fig. [1] indicates the data 



^* ^ vA+A _ vB+B ' W 
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from NSCL/MSU'^. For the experimental free-space NN cross sections, comparing 
theoretical results with experimental data allows us to extract a nuclear symmetry 
energy of i?sym(p) ~ 31.6(p/po)^"^, corresponding to the parameter x = —1 in 
the MDI interaction. With the medium-de penden t NN cross sections, which are 
important for isospin-dependent observables l^" I O" l^ to explain the isospin diffusion 
data requires an ev en s ofter nuclear symmetry energy of -Esym(p) ~ 31.6(p//9o)'' 
with 7 w 0.69 — 1.05!^. In this case, the extracted value for the slope parameter of 
the nuclear symmetry energy at saturation density is L = 88 ±25 MeV and that for 
the isospin-dependent part of the isobaric i ncompres sibility of isospin asymmetric 
nuclear matter is of K^sy — —500 ± 50 MeV 



1451401551 
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Fig. 1. (C olor online) Degree of the isospin diffusio n 1 — Ri as functions of A'asy (left panel, taken 
from Ref. 1^^ and L (right panel, taken from Ref. 1221) using the MDI interaction with x = —2, 
-1, 0, and 1. 



The above extracted symmetry energy agrees with the symmetry energy 



Esymip) — 31.6(p/po) : corresponding to i w 65 MeV and K^, 



-453 MeV, 



that was recently obtained from the isoscaling analyses of isotope ratios in interme- 
diate energy heavy ion collisions'^. The extracted value of i^asy = —500 ±50 MeV 
from the isospin diffusion data is also consistent with the value A'asy = —550 ±100 
MeV obtained from recently measured isotopic dependence of the GMR in even-A 
Sn isotopes '^. These empirically extracted values for L and Ksym represent the 
best current constraints on the nuclear symmetry energy at sub- normal densities. 



4. Constraining effective interactions in the SHF and RMF models 

The non-relativist ic Skyrme-Har tree-Fock (SHF) I 49 | 63 | 64 | 65 | 66 | 67 | g^^^^j relativistic 
mean- field (RMF) I oq I p^ I ' *-* ! ' J- l ' '^ l models are two extensively used phenomcnological 
approaches for nuclear structure studies '■^. In these models, a number of param- 
eters are adjusted to fit the properties of nuclear matter as well as those of many 
finite nuclei. Since there are many different parameter sets in the SHF and RMF 
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models, it is of interest to see to what extend these parameter sets are consistent 
with the constrained symmetry energy from heavy-ion colhsions. 



60 

50 

40 

i. 30 

I20 

UJ 

10 



- Skyrme-Hartree-Fock 
with 21 parameter s 
MDI with x= 
MDi with x^ -1 




0.1 0.2 

p (fm-=) 




0.1 0.2 



Fig. 2. (Color online) Density dependen ce o f Esym{p) for 21 sets of Skyrme interaction 
parameters (left panel, taken from Ref. '^^ and for 23 parameter sets in the nonlinear, 
density-dependent, and point-coupling RMF models (right panel, taken from Ref' — '). 



Fig. [2] displays the density dependence of Esym{p) for 21 sets of Skyrme inter- 
action parameters'^, i.e., SKM, SKM*, RATP, SI, SII, SIII, SIV, SV, SVI, E, E,,, 
G^, R^, Z, Z„, Z;, T, T3, SkX, SkXce, and SkXm (left panel), as well as for 23 
parameter sets in the nonlinear (i.e., NLl^l, ^h2^, NL3'^, NL-SH^Hl, TMl' 
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PKlLLOI^ FSU-GoldL^, HA *"-, NLp l^^LpJ ^ , density-dependent (i.e., TW99 



DD-MEl 1 



DD-ME21 



PKDDl 



DDl 



DD-Fl 



and DDRH-corr 



PC-FII22I, PC-F212»1, PC-F3'221, PC-F4I22I, PC-LAI22I, and 
(righ t panel) . Values of the parameters in the Skyrme 



and poi nt-c oupling (i.e. 

FKVW ESI) RMF models ISni^ ^ 

interactions can be found in Refs. I d-j I p^ I do I pgj. comparison, we also show in Fig. [2] 

results from the MDI interactions with x = —1 (open squares) and (solid squares) 

which give, respectively, the upper and lower boun ds fo r the stiffness of the nuclear 

symmetry energy from the isospin diffusion data 1^^. It is seen that the density 

dependence of the symmetry energy varies drastically among different interactions 

and most of the parameter sets give symmetry energies that are inconsistent w ith 

the extracted one except SIV, SV, Go- , and Ka among the Skyrme parameters '^ 

and TMl, NLp, NLpS, PKDD, and FKVW among the RMF parameter setsl^Sl. 



5. Predictions on the neutron skin thickness of heavy nuclei 

Also affected by the density dependence of nuclear symmetry energy is the neutron 
skin thickness 5 of a nucleus, which is defined as the difference between the root- 
mean-square radii of neutron distribution and proton distribution. In particular, 

the skin thickness of a nucleus is sensitive to the slope parameter L of the nuclear 

■ [7fel5DI51l52l53l45l 



symmetry energy at normal nuclear matter density^ 



Using above 
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21 sets of Skyrme interaction parameters, the neutron skin thickness of several nuclei 
has been studied and it is found that for heavy nuclei there exists a strong hnear 
correlation between S and L. For ^°*Pb. the hnear correlation between S and L is 
given by the following expression 1^^ 

5(2"®Pb) = (0.1066 ± 0.0019) + (0.00133 ± 3.76 x 10"^) x L, (7) 

L = (-78.5 ± 3.2) + (740.4 ± 20.9) x S'(2"®Pb), (8) 

where the units of L and S are MeV and fm. respectively. For ^''^Sn and ^^^Sn, the 
corresponding expressions arel^ 

5("^Sn) = (0.1694 ± 0.0025) + (0.0014 ± 5.12 x 10"^) x L, (9) 

L = (-117.1 ± 5.4) + (695.1 ± 25.3) x SC^^Sn), (10) 

and 

S'(i^''Sn) = (0.1255 ± 0.0020) + (0.0011 ± 4.05 x 10"^) x L, (11) 

L = (-110.1 ± 5.2) + (882.6 ± 32.3) x 5(1243^)^ (-^2) 

Similar linear relations between S and L are also expected for other heavy nuclei. 
This is not surprising since the thickness of the neutron skin in heavy nuclei is 
determined by the pressure difference between n eutron and proton matters in a 
nucleus, which is proportional to the parameter L I < | 49|t)U | oi | o^ | t)3 | 4o | 

The extracted L value from the isospin diffusion data in heavy ion collisions 
allows us to determine from Eqs. (O, (O, and (fTT|) . respectively, a neutron skin 
thickness of 0.22 ± 0.04 fm for ^o^Pb, 0.29 ± 0.04 fm for i^zsn, and 0.22 ± 0.04 fm 
for ^^"^Sn. We note these results are not only in surprisingly good agreement with 
those obtained from the RMF model using an accurately calibrated relativistic 

parametrization, which can describe simultaneously th e g round state properties of 

[701 
finite nuclei and their monopole and dipole resonances Li-^l^ but also consistent with 

the experimental data 1211221931, 

6. Probing the high density behavior of the nuclear symmetry 
energy in heavy-ion reactions 

Although significant progress has been made in the determination of the density 
dependence of the nuclear symmetry energy at sub-normal densities, the high den- 
sity behavior of the the nuclear symmetry energy is still poorly known. Fortunately, 
heavy-ion reactions, especially those induced by high energy radioactive beams to 
be available at high energy radioactive beam facilities, provide a unique opportunity 
to pin down the high density behavior of the symmetry energy. In this section, we 
illustrate via transport model simulations several experimental observables which 
are sensitive to the high density behavior of the symmetry energy. 

Most observables proposed so far for studying the density dependence of the nu- 
clear symmetry energy employ differences or ratios of isospin multiplets of baryons 
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and mes ons as well as mirror nuclei, such as the neutron/proton ratio of emitted 
nuclcons ^^, n eutron -proton differential fl ow l^, neu tron-proton correlat ion tunc- 
tion^, t/3Hel2llMl, ^-/^+ | 22 | 32 | 33 | 36 ] E-/E+l35land K"/K+ ratios IMIMIMI, 

etc. However, some of these observables, especially those involving neutrons, pose 
great challenges in experimental measurements. The measurement of neutrons, par- 
ticularly the low energy ones, always suffers from low detection efficiencies even for 
the most advanced neutron detectors. Therefore, observables involving neutrons 
normally have large systematic errors. Moreover, for essentially all these observ- 
ables, the Coulomb force on charged particles plays an important role and some- 
times competes strongly with the symmetry potential. It is thus very desirable to 
find experimental observables which can reduce the influence of both the Coulomb 
force and the systematic errors associated with the measurement of neutrons. A 
possible candidate for such an observable is the double ratios of emitted particles 
taken from two reaction systems using four isotopes of same clement, such as the 
neutron/proton ratio in the neutron-rich s yste m over that in the more symmetric 
system, as recently done by Famiano et al. '^-'-'. Theoretically, the systematic errors 
associated with transport model calculations, which are mostly related to the uncer- 
tainties in the in-medium NN cross sections, techniques of treating collisions, sizes 
of the lattices in calculating the phase space distributions, techniques in handling 
the Pauli blocking, etc., are also expected to be reduced for relative observables 
that involve double ratios from two similar reaction systems. 
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Fig. 3. (Color online) Double rat io for neutron/proton (left panels, taken from Ref.l^iil) and 7r /7r+ 
(right panel, taken from Ref.ED) in I32gjj_|_l24gjj j^^j Ii2gjj_,_li2gj^ reactions at 400 MeV/nucleon. 



Both the double neutron/proton ratio 1^ and the double 7r^/7r+ ratio 1^ in 
^^^Sn+^^'^Sn and ^^^Sn+^^^Sn reactions at 400 MeV/nucleon have been studied 
recently in the IBUU04 model in order to demonstrate the effect of symmetry 
energy at high density, and the results are shown in Fig. [3l It is seen that these 
ratios have about the same sensitivity to the density dependence of symmetry 
energy as the corresponding single ratio in th e re spective neutron-rich system. Also, 
the double neutron-proton differential flow S2 jn I32gj^_^i24gj^ ^^^^ Ii2gj^_|_ii2gj^ 
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reactions and the n/p ratio of the squeeze-out nuclcons'^in ^^^Sn+^^^'Sn have also 
recently been studied in the IBUU04 model. As shown in Fig. 21 both are indeed 
sensitive to the symmetry energy. In particular, compared to other potential probes, 
the n/p ratio of squeeze-out nucleons carries more direct information about the 
symmetry energy at high densities. The sensitivity to the high density behavior of 
the nuclear symmetry energy observed in the n/p ratio of squeeze-out nucleons with 
high transverse momenta is probably the highest found so far among all observables 
studied within the same transport model. 
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Fig. 4. (Color online) Left panels: Rapidity dependence of the double n-p differential transverse 
flow in reactions of Sn-|-Sn isotopes at 400 and 800 MeV/nucleon with 6 = 5 fm (Taken from 
Ref. Id^. Right panels: Transverse momentum dependence of the ratio of midrapidity neutrons to 
protons emitted in the reaction of ^^^Sn+^'^^Sn at 400 MeV/nucleon and 6 = 5 fm . with (lower 
window) and without (upper window) an azimuthal angle cut (Taken from Ref. l^sl). 



Since the proposal of Aichclin and Ko ^^ that the kaon yield in heavy ion 
collisions at energies that are below the threshold for kaon production in a nucleon- 
nucleon collision in free space may be a sensitive probe of the EOS of nuclear matter 
at high densities, a lot of works have been done both theoretically and experimen- 
tally on this problem '^ ^^ ^ ^-" ^'^ ^ . As the kaon is an iso-doublet meson 
with the quark content of ds for K^ and us for K^ , the K'-* /K^ ratio provides a po- 
tentially good probe of the nuclear symmetry energy as the n/p and tt~ /tt~^ ratios. 
The K^ /K'^ ratio is expected to be particularly sensitive to the high density behav- 
ior of nuclear symmetry energy as kaons are produced mainly from the high density 
region during the early stage of the reaction and thus suffer negligible absorption 
effects. The symmetry energy effect on the K^ /K'^ ratio in heavy-io n colli sions 
has recently been investigated using the UrQMD ^^ and RBUU models I^o I ^p I. The 
results from the RBUU model l^, shown in the left panel of Fig. [5l indicate that 
at beam energies below and around the kinematical threshold of kaon production, 
the K^ / K'^ inclusive yield ratio is more sensitive to the symmetry energy than 
the tt^/tt^. Subthreshold kaon production thus could provide a promising tool to 
extract information on the density dependence of the nuclear symmetry energy. 

Experimentally, the FOPI collaboration at GSI-Darmstadt has recently reported 
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Fig. 5. (Color online) Left panel: it~ /n'^ (upper) and K^ /K^ (lower) ratios as a functi on o f 
the incident energy for central Au+Au collisions with the RBUU model (Taken from Ref. 1^^. 
Right panel: Experimental ratio (K^ / K^)fi^/{K^ / K^)zr (star) and theoretical predictions of 
the thermal model (cross) and the transport model with nucle ar matter (open symbols) and 
realistic (filled symbols) collision scenarios (Taken from Ref. I ^^ j). 



the results on K'^ and K'^ meson production in JfR u 



l^Ru and IIZy 
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collisions at a beam kinetic energy of 1.528 A GeV' I The measured double ratio 
{K^ / K^) Bm / [K^ / K^) Zr is compared in the right panel of Fig.^lto the predictions 
from both the thermal model and the RBUU transport model using two different 
collision scenarios and under different assumptions on the stiffness of the symmetry 
energy. The experimental data show a good agreement with the thermal model 
prediction and also with that of the transport model using the collision scenario 
of an infinite nuclear matter with a soft symmetry energy. Although more realistic 
transport simulations of the collisions show a similar agreement with the data, they 
exhibit a rather weak sensitivity to the symmetry ene rgy. D ue to the complexity 
of subthreshold kaon production in heavy-ion collisions^ ^ I, further experimental 
and theoretical studies are needed to extract useful information on the high density 
behavior of the nuclear symmetry energy from subthreshold kaon production in 
heavy-ion collisions induced by neutron-rich nuclei. 



7. Summary 

Heavy-ion reactions induced by neutron-rich nuclei provide a unique means to in- 
vestigate the equation of state of isospin-asymmctric nuclear matter, especially 
the density dependence of the nuclear symmetry energy. We have reviewed the 
recent progress in extracting the information on the subsaturation density be- 
havior of the nuclear symmetry energy from heavy-ion collisions. From the recent 
analysis of the isospin diffusion data in heavy-ion collisions using an isospin- and 
momentum-dependent transport model with in-medium NN cross sections, a value 
of L = 88 ± 25 MeV for the slope parameter of the nuclear symmetry energy at 
saturation density and a value of ii'asy ~ —500 ± 50 MeV for the isospin-dependent 
part of the isobaric incompressibility of isospin asymmetric nuclear matter have 
been extracted. The extracted symmetry energy agrees with the symmetry energy 
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Esy,n{p) = 31.6(p/po)°■*'^ corresponding to L « 65 MeV and K^^y « -453 MeV, 
that was recently obtained from the isoscaling analyses of isotope ratios in interme- 
diate energy heavy ion collisions "■. The extracted value of /\asy = —500 ±50 MeV 
from the isospin diffusion data is also consistent with the value Kg^sy = —550 ± 100 
MeV obtained from recently measured isotopic dependence of the GMR in even-A 
Sn isotopes '^. These empirically extracted values for L and Ksym represent the 
best current constraints on the nuclear symmetry energy at sub-normal densities. 
This in turn imposes strong constraints on the parameters in both Skyrme and 
RMF effective interactions. Detailed studies have indicated that only SIV, SV, Go-, 
and Rcr among the Skyrme parameters and only TMl, NLp, NLpS, PKDD, and 
FKVW among the RMF parameter sets have symmetry energies that are consis- 
tent with the extracted one. Furthermore, the extracted L value from the isospin 
diffusion data has led to predicted neutron skin thickness of 0.22 ± 0.04 fm for 
208pb, 0.29 ± 0.04 fm for ^^^Sn, and 0.22 ± 0.04 fm for i24Sn. 

We have also reviewed recent theoretical progress in identifying the observables 
in heavy-ion collisions induced by high energy radioactive nuclei that are sensitive 
to the high density behavior of the symmetry energy. Many potentially sensitive 
probes have been found, and they include the tt~ /tt~^ ratio, isospin fractionation, 
n-p differential flow, double n/p and tt~ /n'^ ratios, double n-p differential trans- 
verse flow as well as the K^ /K^ and S~/I]+ ratios. Studying these observable in 
future experiments at high energy radioactive beam facilities is expected to provide 
significant constraints on the behavior of the symmetry energy at supra-normal 
densities. 
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